Transcriptome analysis of the 2,4-dichlorophenoxyacetic acid (2,4-D)-tolerant cotton chromosome substitution line CS-B15sh and its susceptible parental lines G. hirsutum L. cv. Texas Marker-1 and G. barbadense L. cv. Pima 379.

The cotton chromosome substitution line, CS-B15sh, exhibits 41% lower injury from 2,4-D when applied at the field recommended rate of 1.12 kg ae ha-1 (1×) than does Texas Marker-1 (TM-1). CS-B15sh was developed in the genetic background of Gossypium hirsutum L. cv TM-1 and has chromosome introgression on the short arm of chromosome 15 from Gossypium barbadense L. cv. Pima 379. In a previous experiment, we observed reduced translocation of [14C]2,4-D outside the treated leaf tissue in CS-B15sh, which contrasted with an increased translocation of the herbicide in the tissues above and below the treated leaf in TM-1. Our results indicate a potential 2,4-D tolerance mechanism in CS-B15sh involving altered movement of 2,4-D. Here, we used RNA sequencing (RNA-seq) to determine the differential expression of genes between 2,4-D-challenged and control plants of the tolerant (CS-B15sh) and susceptible lines (TM-1 and Pima 379). Several components of the 2,4-D/auxin-response pathway-including ubiquitin E3 ligase, PB1|AUX/IAA, ARF transcription factors, and F-box proteins of the SCFTIR1/AFB complex-were upregulated with at least threefold higher expression in TM-1 compared with CS-B15sh, while both Pima 379 and TM-1 showed the same fold change expression for PB1|AUX/IAA mRNA. Some genes associated with herbicide metabolism, including flavin monooxygenase (Gohir.A01G174100) and FAD-linked oxidase (Gohir.D06G002600), exhibited at least a twofold increase in CS-B15sh than in TM-1 (the gene was not expressed in Pima 379), suggesting a potential relationship between the gene's expression and 2,4-D tolerance. It is interesting to note that glutathione S-transferase was differentially expressed in both CS-B15sh and Pima 379 but not in TM-1, while cytochrome P450 and other genes involved in the oxidation-reduction process were significantly expressed only in CS-B15sh in response to 2,4-D. Gene set enrichment analysis on the union DEGs of the three cotton genotypes revealed the depletion of transcripts involved in photosynthesis and enrichment of transcripts involved in ABA response and signaling.

Exploring the loblolly pine (Pinus taeda L.) genome by BAC sequencing and Cot analysis.

Loblolly pine (LP; Pinus taeda L.) is an economically and ecologically important tree in the southeastern U.S. To advance understanding of the loblolly pine (LP; Pinus taeda L.) genome, we sequenced and analyzed 100 BAC clones and performed a Cot analysis. The Cot analysis indicates that the genome is composed of 57, 24, and 10% highly-repetitive, moderately-repetitive, and single/low-copy sequences, respectively (the remaining 9% of the genome is a combination of fold back and damaged DNA). Although single/low-copy DNA only accounts for 10% of the LP genome, the amount of single/low-copy DNA in LP is still 14 times the size of the Arabidopsis genome. Since gene numbers in LP are similar to those in Arabidopsis, much of the single/low-copy DNA of LP would appear to be composed of DNA that is both gene- and repeat-poor. Macroarrays prepared from a LP bacterial artificial chromosome (BAC) library were hybridized with probes designed from cell wall synthesis/wood development cDNAs, and 50 of the "targeted" clones were selected for further analysis. An additional 25 clones were selected because they contained few repeats, while 25 more clones were selected at random. The 100 BAC clones were Sanger sequenced and assembled. Of the targeted BACs, 80% contained all or part of the cDNA used to target them. One targeted BAC was found to contain fungal DNA and was eliminated from further analysis. Combinations of similarity-based and ab initio gene prediction approaches were utilized to identify and characterize potential coding regions in the 99 BACs containing LP DNA. From this analysis, we identified 154 gene models (GMs) representing both putative protein-coding genes and likely pseudogenes. Ten of the GMs (all of which were specifically targeted) had enough support to be classified as intact genes. Interestingly, the 154 GMs had statistically indistinguishable (α = 0.05) distributions in the targeted and random BAC clones (15.18 and 12.61 GM/Mb, respectively), whereas the low-repeat BACs contained significantly fewer GMs (7.08 GM/Mb). However, when GM length was considered, the targeted BACs had a significantly greater percentage of their length in GMs (3.26%) when compared to random (1.63%) and low-repeat (0.62%) BACs. The results of our study provide insight into LP evolution and inform ongoing efforts to produce a reference genome sequence for LP, while characterization of genes involved in cell wall production highlights carbon metabolism pathways that can be leveraged for increasing wood production.

Proteomic analysis of the venom and venom sac of the woodwasp, Sirex noctilio - Towards understanding its biological impact.

UNLABELLED: The European horntail woodwasp, Sirex noctilio, is an invasive insect that attacks conifer hosts, particularly Pinus species. Venom injected by female S. noctilio, together with its symbiotic fungus, damages the normal physiology of Pinus, leading to death of the tree. To identify the proteinaceous components in the venom and uncover the interplay between venom proteins and tree proteins, clarification of the overall profile of proteins produced in the venom gland apparatus was carried out in this work. The venom sac proteome utilised in-solution digested in either a natural or deglycosylated state, prior to nanoHPLC LTQ-Orbitrap under CID/ETD mode. Here, we report the identification of 1454 and 1225 proteins in venom and sac, respectively, with 410 mutual proteins. Approximately 90 proteins were predicted to be secretory, of which 8 have features characteristic of toxins. Chemosensory binding proteins were also identified. Gene ontology and KEGG pathway analysis were employed to predict the protein functions and biological pathways in venom and sac. Protein-protein interaction (PPI) analysis suggested that one-step responses represent the majority of the Sirex-Pinus PPIs, and the proteins representing network hub nodes could be of importance for the development of pest management strategies. SIGNIFICANCE: The woodwasp Sirex noctilio is an invasive species in many parts of the world, including Australia and North America, where it is considered within the top 10 most serious forest insects. Where they have been introduced, the female woodwasps attack living pine trees, causing significant economic losses. Central to this destruction is the woodwasp's life cycle requirement to bore a hole to deposit eggs and a toxic mucus that disables the tree's network for transporting water and nutrients, yet aids in larval survival. Here we specifically examine the mucus gland apparatus and its contents, revealing the protein components that together with 'noctilisin' facilitate this complex association. The identification of chemosensory binding proteins further supports a role for the woodwasp ovipositor as an instrument for early stages of host tree selection. These findings could provide important clues towards the development of novel control tools against this pest.

Activation of defence pathways in Scots pine bark after feeding by pine weevil (Hylobius abietis).

BACKGROUND: During their lifetime, conifer trees are exposed to numerous herbivorous insects. To protect themselves against pests, trees have developed a broad repertoire of protective mechanisms. Many of the plant's defence reactions are activated upon an insect attack, and the underlying regulatory mechanisms are not entirely understood yet, in particular in conifer trees. Here, we present the results of our studies on the transcriptional response and the volatile compounds production of Scots pine (Pinus sylvestris) upon the large pine weevil (Hylobius abietis) feeding. RESULTS: Transcriptional response of Scots pine to the weevil attack was investigated using a novel customised 36.4 K Pinus taeda microarray. The weevil feeding caused large-scale changes in the pine transcriptome. In total, 774 genes were significantly up-regulated more than 4-fold (p≤0.05), whereas 64 genes were significantly down-regulated more than 4-fold. Among the up-regulated genes, we could identify genes involved in signal perception, signalling pathways, transcriptional regulation, plant hormone homeostasis, secondary metabolism and defence responses. The weevil feeding on stem bark of pine significantly increased the total emission of volatile organic compounds from the undamaged stem bark area. The emission levels of monoterpenes and sesquiterpenes were also increased. Interestingly, we could not observe any correlation between the increased production of the terpenoid compounds and expression levels of the terpene synthase-encoding genes. CONCLUSIONS: The obtained data provide an important insight into the transcriptional response of conifer trees to insect herbivory and illustrate the massive changes in the host transcriptome upon insect attacks. Moreover, many of the induced pathways are common between conifers and angiosperms. The presented results are the first ones obtained by the use of a microarray platform with an extended coverage of pine transcriptome (36.4 K cDNA elements). The platform will further facilitate the identification of resistance markers with the direct relevance for conifer tree breeding.

Decoding the massive genome of loblolly pine using haploid DNA and novel assembly strategies.

BACKGROUND: The size and complexity of conifer genomes has, until now, prevented full genome sequencing and assembly. The large research community and economic importance of loblolly pine, Pinus taeda L., made it an early candidate for reference sequence determination. RESULTS: We develop a novel strategy to sequence the genome of loblolly pine that combines unique aspects of pine reproductive biology and genome assembly methodology. We use a whole genome shotgun approach relying primarily on next generation sequence generated from a single haploid seed megagametophyte from a loblolly pine tree, 20-1010, that has been used in industrial forest tree breeding. The resulting sequence and assembly was used to generate a draft genome spanning 23.2 Gbp and containing 20.1 Gbp with an N50 scaffold size of 66.9 kbp, making it a significant improvement over available conifer genomes. The long scaffold lengths allow the annotation of 50,172 gene models with intron lengths averaging over 2.7 kbp and sometimes exceeding 100 kbp in length. Analysis of orthologous gene sets identifies gene families that may be unique to conifers. We further characterize and expand the existing repeat library based on the de novo analysis of the repetitive content, estimated to encompass 82% of the genome. CONCLUSIONS: In addition to its value as a resource for researchers and breeders, the loblolly pine genome sequence and assembly reported here demonstrates a novel approach to sequencing the large and complex genomes of this important group of plants that can now be widely applied.

Noctilisin, a Venom Glycopeptide of Sirex noctilio (Hymenoptera: Siricidae), Causes Needle Wilt and Defense Gene Responses in Pines.

During oviposition, female Sirex noctilio (F.) (Siricidae) woodwasps inject their conifer hosts with a venom gland secretion. The secretion induces a variety of host physiological changes that facilitate subsequent lethal infection by a symbiotic fungus. A heat-stable factor that can migrate from the site of oviposition in the trunk through the xylem to needles in the crown of attacked pines was purified by size-fractionation and reversed-phase-high-performance liquid chromatography using activity assays based on defense gene induction as well as the needle wilt response in pine shoot explants. An 11-amino acid, posttranslationally modified peptide (SEGPROGTKRP) encoded by the most abundant transcript recovered from S. noctilio venom gland tissue comprised the backbone of the 1,850 Da active factor. Posttranslational modifications included hydroxylation of a Pro residue at position 6 as well as O-glycosylation of Ser and Thr residues at positions 1 and 8, respectively. The O-linked sugars were identical α-linked N-acetylgalactosamine residues modified at the C6 position by addition of phosphoethanolamine. In contrast to the native peptide, a synthetic version of the hydroxylated peptide backbone lacking the glycosyl side chains failed to induce pine defense genes or cause needle wilt in excised shoots. This peptide, hereafter called noctilisin, is related to the O-glycosylated short-chain proline-rich antimicrobial peptides exemplified by drosocin. The noctilisin structure contains motifs which may explain how it avoids detection by pine defense systems.

Transcriptomic analysis highlights epigenetic and transcriptional regulation during zygotic embryo development of Pinus pinaster.

BACKGROUND: It is during embryogenesis that the plant body plan is established and the meristems responsible for all post-embryonic growth are specified. The molecular mechanisms governing conifer embryogenesis are still largely unknown. Their elucidation may contribute valuable information to clarify if the distinct features of embryo development in angiosperms and gymnosperms result from differential gene regulation. To address this issue, we have performed the first transcriptomic analysis of zygotic embryo development in a conifer species (Pinus pinaster) focusing our study in particular on regulatory genes playing important roles during plant embryo development, namely epigenetic regulators and transcription factors. RESULTS: Microarray analysis of P. pinaster zygotic embryogenesis was performed at five periods of embryo development from early developing to mature embryos. Our results show that most changes in transcript levels occurred in the first and the last embryo stage-to-stage transitions, namely early to pre-cotyledonary embryo and cotyledonary to mature embryo. An analysis of functional categories for genes that were differentially expressed through embryogenesis highlighted several epigenetic regulation mechanisms. While putative orthologs of transcripts associated with mechanisms that target transposable elements and repetitive sequences were strongly expressed in early embryogenesis, PRC2-mediated repression of genes seemed more relevant during late embryogenesis. On the other hand, functions related to sRNA pathways appeared differentially regulated across all stages of embryo development with a prevalence of miRNA functions in mid to late embryogenesis. Identification of putative transcription factor genes differentially regulated between consecutive embryo stages was strongly suggestive of the relevance of auxin responses and regulation of auxin carriers during early embryogenesis. Such responses could be involved in establishing embryo patterning. Later in development, transcripts with homology to genes acting on modulation of auxin flow and determination of adaxial-abaxial polarity were up-regulated, as were putative orthologs of genes required for meristem formation and function as well as establishment of organ boundaries. Comparative analysis with A. thaliana embryogenesis also highlighted genes involved in auxin-mediated responses, as well as epigenetic regulation, indicating highly correlated transcript profiles between the two species. CONCLUSIONS: This is the first report of a time-course transcriptomic analysis of zygotic embryogenesis in a conifer. Taken together our results show that epigenetic regulation and transcriptional control related to auxin transport and response are critical during early to mid stages of pine embryogenesis and that important events during embryogenesis seem to be coordinated by putative orthologs of major developmental regulators in angiosperms.

A SNP resource for Douglas-fir: de novo transcriptome assembly and SNP detection and validation.

BACKGROUND: Douglas-fir (Pseudotsuga menziesii), one of the most economically and ecologically important tree species in the world, also has one of the largest tree breeding programs. Although the coastal and interior varieties of Douglas-fir (vars. menziesii and glauca) are native to North America, the coastal variety is also widely planted for timber production in Europe, New Zealand, Australia, and Chile. Our main goal was to develop a SNP resource large enough to facilitate genomic selection in Douglas-fir breeding programs. To accomplish this, we developed a 454-based reference transcriptome for coastal Douglas-fir, annotated and evaluated the quality of the reference, identified putative SNPs, and then validated a sample of those SNPs using the Illumina Infinium genotyping platform. RESULTS: We assembled a reference transcriptome consisting of 25,002 isogroups (unique gene models) and 102,623 singletons from 2.76 million 454 and Sanger cDNA sequences from coastal Douglas-fir. We identified 278,979 unique SNPs by mapping the 454 and Sanger sequences to the reference, and by mapping four datasets of Illumina cDNA sequences from multiple seed sources, genotypes, and tissues. The Illumina datasets represented coastal Douglas-fir (64.00 and 13.41 million reads), interior Douglas-fir (80.45 million reads), and a Yakima population similar to interior Douglas-fir (8.99 million reads). We assayed 8067 SNPs on 260 trees using an Illumina Infinium SNP genotyping array. Of these SNPs, 5847 (72.5%) were called successfully and were polymorphic. CONCLUSIONS: Based on our validation efficiency, our SNP database may contain as many as ~200,000 true SNPs, and as many as ~69,000 SNPs that could be genotyped at ~20,000 gene loci using an Infinium II array-more SNPs than are needed to use genomic selection in tree breeding programs. Ultimately, these genomic resources will enhance Douglas-fir breeding and allow us to better understand landscape-scale patterns of genetic variation and potential responses to climate change.

Biomarker genes highlight intraspecific and interspecific variations in the responses of Pinus taeda L. and Pinus radiata D. Don to Sirex noctilio F. acid gland secretions.

Sirex noctilio F., a Eurasian horntail woodwasp recently introduced into North America, oviposits in pines and other conifers and in the process spreads a phytopathogenic fungus that serves as a food source for its larvae. During oviposition the woodwasp also deposits mucus produced in its acid (venom) gland that alters pine defense responses and facilitates infection by the fungus. A 26,496-feature loblolly pine cDNA microarray was used to survey gene expression of pine tissue responding to S. noctilio venom. Six genes were selected for further assessment by quantitative real-time polymerase chain reaction (qRT-PCR), including one that encoded an apparent PR-4 protein and another that encoded a thaumatin-like protein. Expression of both was strongly induced in response to venom, while expression of an apparent actin gene (ACT1) was stable in response to the venom. The pattern of gene response was similar in Pinus taeda L. and Pinus radiata D. Don, but the magnitude of response in P. radiata was significantly stronger for each of the induced genes. The magnitude of the biomarker gene response to venom also varied according to genotype within these two species. The qRT-PCR assay was used to demonstrate that the primary bioactive component in S. noctilio venom is a polypeptide.

Towards decoding the conifer giga-genome.

Several new initiatives have been launched recently to sequence conifer genomes including pines, spruces and Douglas-fir. Owing to the very large genome sizes ranging from 18 to 35 gigabases, sequencing even a single conifer genome had been considered unattainable until the recent throughput increases and cost reductions afforded by next generation sequencers. The purpose of this review is to describe the context for these new initiatives. A knowledge foundation has been acquired in several conifers of commercial and ecological interest through large-scale cDNA analyses, construction of genetic maps and gene mapping studies aiming to link phenotype and genotype. Exploratory sequencing in pines and spruces have pointed out some of the unique properties of these giga-genomes and suggested strategies that may be needed to extract value from their sequencing. The hope is that recent and pending developments in sequencing technology will contribute to rapidly filling the knowledge vacuum surrounding their structure, contents and evolution. Researchers are also making plans to use comparative analyses that will help to turn the data into a valuable resource for enhancing and protecting the world's conifer forests.

The phenylalanine ammonia lyase (PAL) gene family shows a gymnosperm-specific lineage.

BACKGROUND: Phenylalanine ammonia lyase (PAL) is a key enzyme of the phenylpropanoid pathway that catalyzes the deamination of phenylalanine to trans-cinnamic acid, a precursor for the lignin and flavonoid biosynthetic pathways. To date, PAL genes have been less extensively studied in gymnosperms than in angiosperms. Our interest in PAL genes stems from their potential role in the defense responses of Pinus taeda, especially with respect to lignification and production of low molecular weight phenolic compounds under various biotic and abiotic stimuli. In contrast to all angiosperms for which reference genome sequences are available, P. taeda has previously been characterized as having only a single PAL gene. Our objective was to re-evaluate this finding, assess the evolutionary history of PAL genes across major angiosperm and gymnosperm lineages, and characterize PAL gene expression patterns in Pinus taeda. METHODS: We compiled a large set of PAL genes from the largest transcript dataset available for P. taeda and other conifers. The transcript assemblies for P. taeda were validated through sequencing of PCR products amplified using gene-specific primers based on the putative PAL gene assemblies. Verified PAL gene sequences were aligned and a gene tree was estimated. The resulting gene tree was reconciled with a known species tree and the time points for gene duplication events were inferred relative to the divergence of major plant lineages. RESULTS: In contrast to angiosperms, gymnosperms have retained a diverse set of PAL genes distributed among three major clades that arose from gene duplication events predating the divergence of these two seed plant lineages. Whereas multiple PAL genes have been identified in sequenced angiosperm genomes, all characterized angiosperm PAL genes form a single clade in the gene PAL tree, suggesting they are derived from a single gene in an ancestral angiosperm genome. The five distinct PAL genes detected and verified in P. taeda were derived from a combination of duplication events predating and postdating the divergence of angiosperms and gymnosperms. CONCLUSIONS: Gymnosperms have a more phylogenetically diverse set of PAL genes than angiosperms. This inference has contrasting implications for the evolution of PAL gene function in gymnosperms and angiosperms.

Microarray analysis and scale-free gene networks identify candidate regulators in drought-stressed roots of loblolly pine (P. taeda L.).

BACKGROUND: Global transcriptional analysis of loblolly pine (Pinus taeda L.) is challenging due to limited molecular tools. PtGen2, a 26,496 feature cDNA microarray, was fabricated and used to assess drought-induced gene expression in loblolly pine propagule roots. Statistical analysis of differential expression and weighted gene correlation network analysis were used to identify drought-responsive genes and further characterize the molecular basis of drought tolerance in loblolly pine. RESULTS: Microarrays were used to interrogate root cDNA populations obtained from 12 genotype × treatment combinations (four genotypes, three watering regimes). Comparison of drought-stressed roots with roots from the control treatment identified 2445 genes displaying at least a 1.5-fold expression difference (false discovery rate = 0.01). Genes commonly associated with drought response in pine and other plant species, as well as a number of abiotic and biotic stress-related genes, were up-regulated in drought-stressed roots. Only 76 genes were identified as differentially expressed in drought-recovered roots, indicating that the transcript population can return to the pre-drought state within 48 hours. Gene correlation analysis predicts a scale-free network topology and identifies eleven co-expression modules that ranged in size from 34 to 938 members. Network topological parameters identified a number of central nodes (hubs) including those with significant homology (E-values ≤ 2 × 10-30) to 9-cis-epoxycarotenoid dioxygenase, zeatin O-glucosyltransferase, and ABA-responsive protein. Identified hubs also include genes that have been associated previously with osmotic stress, phytohormones, enzymes that detoxify reactive oxygen species, and several genes of unknown function. CONCLUSION: PtGen2 was used to evaluate transcriptome responses in loblolly pine and was leveraged to identify 2445 differentially expressed genes responding to severe drought stress in roots. Many of the genes identified are known to be up-regulated in response to osmotic stress in pine and other plant species and encode proteins involved in both signal transduction and stress tolerance. Gene expression levels returned to control values within a 48-hour recovery period in all but 76 transcripts. Correlation network analysis indicates a scale-free network topology for the pine root transcriptome and identifies central nodes that may serve as drivers of drought-responsive transcriptome dynamics in the roots of loblolly pine.

Differential responses of the promoters from nearly identical paralogs of loblolly pine (Pinus taeda L.) ACC oxidase to biotic and abiotic stresses in transgenic Arabidopsis thaliana.

Promoters from an ACC oxidase gene (PtACO1) and its nearly identical paralog (NIP) (PtACO2) of loblolly pine (Pinus taeda L.) were recovered from genomic DNA using PCR amplification. Transgenic Arabidopsis plants harboring genetic constructs from which beta-glucuronidase (GUS) expression was driven by the full-length (pACO1:GUS, pACO2:GUS) or truncated (pACO1-1.2:GUS, pACO2-1.2:GUS) loblolly pine ACC oxidase gene promoters displayed distinctive patterns of expression for the different promoter constructs. Both full-length promoter constructs, but not those using truncated promoters, responded to indole-3-acetic acid (IAA) and wounding. Both pACO1:GUS and pACO1-1.2:GUS responded to pathogen attack, while neither version of the pACO2 promoter responded to infection. In the inflorescence stalks, the full-length pACO1 promoter construct, but not the truncated pACO1-1.2:GUS or either pACO2 construct, responded to bending stress. When flowering transgenic Arabidopsis plants were placed in a horizontal position for 48 h, expression from pACO2:GUS, but not the other constructs, was induced on the underside of shoots undergoing gravitropic reorientation. The expression pattern for the pACO2:GUS construct in transgenic Arabidopsis was consistent with what might be expected for a gene promoter involved in the compression wood response in loblolly pine. Although near complete sequence identity between PtACO1 and PtACO2 transcripts prevented quantitation of specific gene products, the promoter expression analyses presented in this study provide strong evidence that the two ACC oxidase genes are likely differentially expressed and responded to different external stimuli in pine. These results are discussed with respect to the potential functional differences between these two genes in loblolly pine.

An improved method of RNA isolation from loblolly pine (P. taeda L.) and other conifer species.

Tissues isolated from conifer species, particularly those belonging to the Pinaceae family, such as loblolly pine (Pinus taeda L.), contain high concentrations of phenolic compounds and polysaccharides that interfere with RNA purification. Isolation of high-quality RNA from these species requires rigorous tissue collection procedures in the field and the employment of an RNA isolation protocol comprised of multiple organic extraction steps in order to isolate RNA of sufficient quality for microarray and other genomic analyses. The isolation of high-quality RNA from field-collected loblolly pine samples can be challenging, but several modifications to standard tissue and RNA isolation procedures greatly improve results. The extent of general RNA degradation increases if samples are not properly collected and transported from the field, especially during large-scale harvests. Total RNA yields can be increased significantly by pulverizing samples in a liquid nitrogen freezer mill prior to RNA isolation, especially when samples come from woody tissues. This is primarily due to the presence of oxidizing agents, such as phenolic compounds, and polysaccharides that are both present at high levels in extracts from the woody tissues of most conifer species. If not removed, these contaminants can carry over leading to problems, such as RNA degradation, that result in low yields and a poor quality RNA sample. Carryover of phenolic compounds, as well as polysaccharides, can also reduce or even completely eliminate the activity of reverse transcriptase or other polymerases commonly used for cDNA synthesis. In particular, RNA destined to be used as template for double-stranded cDNA synthesis in the generation of cDNA libraries, single-stranded cDNA synthesis for PCR or qPCR's, or for the synthesis of microarray target materials must be of the highest quality if researchers expect to obtain optimal results. RNA isolation techniques commonly employed for many other plant species are often insufficient in their ability to remove these contaminants from conifer samples and thus do not yield total RNA samples suitable for downstream manipulations. In this video we demonstrate methods for field collection of conifer tissues, beginning with the felling of a forty year-old tree, to the harvesting of phloem, secondary xylem, and reaction wood xylem. We also demonstrate an RNA isolation protocol that has consistently yielded high-quality RNA for subsequent enzymatic manipulations.

Processing the loblolly pine PtGen2 cDNA microarray.

PtGen2 is a 26,496 feature cDNA microarray containing amplified loblolly pine ESTs. The array is produced in our laboratory for use by researchers studying gene expression in pine and other conifer species. PtGen2 was developed as a result of our gene discovery efforts in loblolly pine, and is comprised of sequences identified primarily from root tissues, but also from needle and stem. PtGen2 has been tested by hybridizing different Cy-dye labeled conifer target cDNAs, using both amplified and non-amplified indirect labeling methods, and also tested with a number of hybridization and washing conditions. This video focuses on the handling and processing of slides before and after pre-hybridization, as well as after hybridization, using some modifications to procedures developed previously. Also included, in text form only, are the protocols used for the generation, labeling and clean up of target cDNA s, as well as information on software used for downstream data processing. PtGen2 is printed with a proprietary print buffer that contains high concentrations of salt that can be difficult to remove completely. The slides are washed first in a warm SDS solution prior to pre-hybridization. After pre-hybridization, the slides are washed vigorously in several changes of water to complete removal of remaining salts. LifterSlips are then cleaned and positioned on the slides and labeled cDNA is carefully loaded onto the microarray by way of capillary action which provides for even distribution of the sample across the slide, and reduces the chance of bubble incorporation. Hybridization of targets to the array is done at 48 degrees C in high humidity conditions. After hybridization, a series of standard washes are done at 53 degrees C and room temperature for extended times. Processing PtGen2 slides using this technique reduces salt and SDS-derived artifacts often seen when the array is processed less rigorously. Hybridizing targets derived from several different conifer RNA sources, this processing protocol yielded fewer artifacts, reduced background, and provided better consistency among different experimental groups of arrays.

Carbohydrate-related genes and cell wall biosynthesis in vascular tissues of loblolly pine (Pinus taeda).

Loblolly pine (Pinus taeda L.), the most widely planted tree species in the United States, is an important source of wood and wood fibers for a multitude of consumer products. Wood fibers are primarily composed of secondary cell walls, and cellulose, hemicelluloses and lignin are major components of wood. Fiber morphology and cell wall composition are important determinants of wood properties. We used comparative genomics to identify putative genes for cellulose and hemicellulose synthesis in loblolly pine that are homologous to genes implicated in cell wall synthesis in angiosperms. Sequences encoding putative secondary cell wall cellulose synthase genes, cellulose synthase-like genes, a membrane-bound endoglucanase gene, a sucrose synthase gene, a UDP-glucose pyrophosphorylase gene and GDP-mannose pyrophosphorylase genes were identified in expressed sequence tag (EST) collections from loblolly pine. Full-length coding sequences were obtained from cDNA clones isolated from a library constructed from developing xylem. Phylogenetic relationships between the genes from loblolly pine and angiosperm taxa were examined and transcriptional profiling in vascular tissues was conducted by real-time quantitative, reverse transcriptase-polymerase chain reaction. The putative cell wall synthesis genes were expressed at high levels in vascular tissues and a subset was differentially regulated in xylem and phloem tissues. Inferred phylogenetic relationships and expression patterns for the genes from loblolly pine were consistent with roles in synthesis of complex carbohydrates of the cell wall. These studies suggest functional conservation of homologous wood formation genes in gymnosperm and angiosperm taxa.

WebTraceMiner: a web service for processing and mining EST sequence trace files.

Expressed sequence tags (ESTs) remain a dominant approach for characterizing the protein-encoding portions of various genomes. Due to inherent deficiencies, they also present serious challenges for data quality control. Before GenBank submission, EST sequences are typically screened and trimmed of vector and adapter/linker sequences, as well as polyA/T tails. Removal of these sequences presents an obstacle for data validation of error-prone ESTs and impedes data mining of certain functional motifs, whose detection relies on accurate annotation of positional information for polyA tails added posttranscriptionally. As raw DNA sequence information is made increasingly available from public repositories, such as NCBI Trace Archive, new tools will be necessary to reanalyze and mine this data for new information. WebTraceMiner (www.conifergdb.org/software/wtm) was designed as a public sequence processing service for raw EST traces, with a focus on detection and mining of sequence features that help characterize 3' and 5' termini of cDNA inserts, including vector fragments, adapter/linker sequences, insert-flanking restriction endonuclease recognition sites and polyA or polyT tails. WebTraceMiner complements other public EST resources and should prove to be a unique tool to facilitate data validation and mining of error-prone ESTs (e.g. discovery of new functional motifs).

ConiferEST: an integrated bioinformatics system for data reprocessing and mining of conifer expressed sequence tags (ESTs).

BACKGROUND: With the advent of low-cost, high-throughput sequencing, the amount of public domain Expressed Sequence Tag (EST) sequence data available for both model and non-model organism is growing exponentially. While these data are widely used for characterizing various genomes, they also present a serious challenge for data quality control and validation due to their inherent deficiencies, particularly for species without genome sequences. DESCRIPTION: ConiferEST is an integrated system for data reprocessing, visualization and mining of conifer ESTs. In its current release, Build 1.0, it houses 172,229 loblolly pine EST sequence reads, which were obtained from reprocessing raw DNA sequencer traces using our software--WebTraceMiner. The trace files were downloaded from NCBI Trace Archive. ConiferEST provides biologists unique, easy-to-use data visualization and mining tools for a variety of putative sequence features including cloning vector segments, adapter sequences, restriction endonuclease recognition sites, polyA and polyT runs, and their corresponding Phred quality values. Based on these putative features, verified sequence features such as 3' and/or 5' termini of cDNA inserts in either sense or non-sense strand have been identified in-silico. Interestingly, only 30.03% of the designated 3' ESTs were found to have an authenticated 5' terminus in the non-sense strand (i.e., polyT tails), while fewer than 5.34% of the designated 5' ESTs had a verified 5' terminus in the sense strand. Such previously ignored features provide valuable insight for data quality control and validation of error-prone ESTs, as well as the ability to identify novel functional motifs embedded in large EST datasets. We found that "double-termini adapters" were effective indicators of potential EST chimeras. For all sequences with in-silico verified termini/terminus, we used InterProScan to assign protein domain signatures, results of which are available for in-depth exploration using our biologist-friendly web interfaces. CONCLUSION: ConiferEST represents a unique and complementary public resource for EST data integration and mining in conifers by reprocessing raw DNA traces, identifying putative sequence features and determining and annotating in-silico verified features. Seamlessly integrated with other public resources, ConiferEST provides biologists powerful tools to verify data, visualize abnormalities, including EST chimeras, and explore large EST datasets.

Water stress-responsive genes in loblolly pine (Pinus taeda) roots identified by analyses of expressed sequence tag libraries.

Drought stress is the principal cause of seedling mortality in pine forests of the southeastern United States and in many other forested regions around the globe. As part of a larger effort to discover loblolly pine genes, this study subjected rooted cuttings of three unrelated pine genotypes to three watering regimens. Expressed sequence tags (ESTs) were obtained from both the 3' and 5' ends of 12,918 randomly selected cDNAs generated from root tissues. These ESTs were clustered to identify 6,765 unique transcripts (UniScripts) derived from 6,202 putative unique genes (UniGenes-S). Tentative annotations were assigned on the basis of BLASTX comparisons to the Protein Information Resource Nonredundant Reference (PIR-NREF) database. Expression levels of 42 UniScripts varied with high statistical significance with respect to treatment. Many of them resembled gene products shown to be important for drought tolerance in other species, including dehydrins, endochitinases, cytochrome P450 enzymes, pathogenesis-related proteins and various late-embryogenesis abundant (LEA) gene products. Similarly, expression levels of 110 UniScripts varied with high statistical significance among genotypes, indicating that gene expression patterns in this species are much more dependent on genotype than on treatment. Most of the water stress-induced pine UniScripts that appeared to encode products resembling drought tolerance factors in other species were most highly induced in a single genotype, suggesting that particularly useful adaptive alleles for drought tolerance might exist within the collection of cDNAs characterized from this genotype. Mining and visualizing the complete data set, as well as downloading of both EST and UniScript contig sequences, are possible using MAGIC Gene Discovery at http://fungen.org/genediscovery/.

Gene expression changes in Arabidopsis thaliana seedling roots exposed to the munition hexahydro-1,3,5-trinitro-1,3,5-triazine.

Arabidopsis thaliana root transcriptome responses to the munition, hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX), were assessed using serial analysis of gene expression (SAGE). Sequencing of SAGE libraries from control and RDX-exposed root tissues revealed induction of genes known to respond to a variety of general stresses. Among the highly induced genes were several encoding molecular chaperones and transcription factors as well as vacuolar proteins and peroxidases. Strongly repressed transcripts included ones encoding ribosomal proteins, a cyclophilin, a katanin, and a peroxidase. Comparison of the transcriptional profile for the RDX response to a profile previously described for Arabidopsis roots exposed to trinitrotoluene (TNT) revealed significant differences in the inferred gene expression patterns. This suggests that Arabidopsis employs drastically different mechanisms for coping with these two compounds. With respect to the goal of engineering plants to better tolerate and degrade explosives at contaminated sites, these results suggest that enhancement of different genes and metabolic pathways may be required to deal effectively with each type of explosive. This has ramifications for phytoremediation efforts since many contaminated sites harbor both compounds.